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We propose a general formalism for analytical description of multiatomic ensembles interacting 
with a single mode quantized cavity field under the assumption that most atoms remain un-excited 
on average. By combining the obtained formalism with the nilpotent technique for the description of 
multipartite entanglement we are able to overview in a unified fashion different probabilistic control 
scenarios of entanglement among atoms or examine atomic ensembles. We then apply the proposed 
control schemes to the creation of multiatom states useful for quantum information. 
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I. INTRODUCTION 

Engineering of the entanglement properties of multi- 
partite states has recently become the subject of exten- 
sive research due to its relevance to quantum information 
processing and computation. A number of entanglement 
control schemes based on neutral atoms in a quantized 
cavity field [l|, 0, H, 0] or in optical lattices [Sj] as well 
as trapped-ion setups Q have been proposed, and im- 
plemented experimentally. However, a universal theory 
of multipartite entanglement, which could guide the di- 
verse experimental and theoretical efforts, is still missing. 
Instead, most of these proposals are focused on the gen- 
eration of specific entangled states rather than general 
states with an arbitrary chosen entanglement. 

Here, we seek a general framework for multipartite en- 
tangled state engineering of neutral atoms in a quantized 
cavity by making use of the recently introduced formal- 
ism for the exhaustive description multipartite entangle- 
ment based on the nilpotent polynomial technique [7]. 

The nilpotent formalism relies on raising operators act- 
ing on a reference product state-for two-level systems 
the variables of these polynomials are the nilpotent op- 
erators . Every quantum state can be represented 
as a polynomial of operators F[{(7^}) acting on a ref- 
erence vector which we can choose to be the "atomic 
vacuum" \0) = ITl'^)?!- main object of interest 

n 

is the logarithm of F{{5f}), that we call "nilpotential" 
/ = h\F{{a^}). It provides us with a simple criterion of 
entanglement for any binary partition, A and B, of the 
quantum system: the subsystems A and B arc disentan- 
gled iff Jaub = fA + Ib- By the term "tanglemeter" /c 
we denote the nilpotential of the "canonic state" which is 
the closest to the atomic vacuum state among the variety 
of the states subject to all possible local transformations. 
This polynomial contains exhaustive information about 
the entanglement. 

In order to implement the nilpotent polynomial tech- 
nique to describe controlled entanglement among TV neu- 
tral two-level atoms interacting with a quantized single- 



mode electromagnetic field in a high-Q cavity, we further 
assume that the number of excited atoms remains low 
on average, well below N, during the interaction. In 
other words, we consider an entangled state of a weakly 
excited multiatom ensemble that is coupled to the quan- 
tized field. This assumption allows us to obtain an exact 
analytic description, which can be directly interpreted in 
terms of the nilpotent formalism for multipartite entan- 
glement, and propose several techniques for rather gen- 
eral multiatom entanglement engineering. 

Unitary control over the system is exerted by squeez- 
ing, displacing, or Kerr-nonlinear self-modulation of the 
cavity field, as well as manipulating the coupling between 
the atoms and the field moving the atoms in the cavity. 
This unitary control, combined with nonunitary measure- 
ments of the number of cavity photons, that is projec- 
tion of the cavity-field state onto a photon number state, 
yields a vast variety of possible outcomes. We detail such 
a probabilistic approach to the construction of the Dicke 
state for N two-level atoms with M excitations p, 's'l , the 
W [H and GHZ [13] states. We also present a proce- 
dure for entangling two atomic ensembles. Decoherence 
is ignored throughout the analysis. This is justified by 
the fact that in our analysis high-Q cavities are required 
and also that the atoms remain low-excited on average. 

The structure of the paper is the following. In SecHIlwe 
derive the analytic description of weakly excited atomic 
ensembles. In Sec. IIIII we turn to the entanglement con- 
trol schemes, which we further specify upon discussing 
their particular applications in Sec. [TV] We conclude by 
discussing the results obtained. Details of the calcula- 
tions are included in Appendices A and B. 



II. ANALYTIC DESCRIPTION OF 
TWO-LEVELS ATOMS IN A CAVITY 

In this section we derive an analytic expression for the 
time-dependent entangled state of N atoms and the cav- 
ity photons,with the help of the functional integration 
technique, assuming that the number of atoms in the ex- 
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cited state remains small as compared to A''. 

The Hamiltonian of a single-mode cavity field coupled 
to N two-level atoms via laser-induced Raman interac- 
tion reads, in the interaction representation, 

n 

(1) 

Here ^md w„ are the resonant frequencies of the cavity 
mode and the n-th atom, respectively, C„(t) is the con- 
trolled Raman coupling of the n-th atom with the cavity 
field induced by an external laser field £{t), ct^'"^'" are 
the Pauli matrices, and a\ a are the field-mode creation 
and annihilation operators, respectively. 

On separating the atomic and cavity field variables 
with the help of a functional integral over the complex 
functional variable z{t), we express the evolution opera- 
tor as 



Vz{t)Vz*it) 
U = / exp 



-i \z{t)rdt 



N 



Uoit) n Unit), 

(2) 



n=l 



where the evolution operators Un,o{i) in the functional 
integral satisfy the dynamic equations 



i^C/oW = [Lom\a-E{t)z*{t){a\+a)]Ua(t) (3) 

i|f7„(i) = [^d^^ + z{t)a,{t) (a+ + d-)] Ur.it), (4) 



and 



^ = J Vz{t)Vz*{t)exp -i J \z{t)\'^dt 



(5) 



is a normalization factor. 

For the field and the atoms initially in the ground 
states (|0)p and \0)^, respectively), the Schrodinger equa- 
tion gives 



IV'o (t)) = exp 



t t t 

iat y"e''^«(^"*)£-(T)z*(T)dr-y" j e''^^'''^-''^ E{t)E{9)z* {t)z* {e)dTde 

T 



|0)„ , 



(6) 



while the second-order time-dependent perturbation theory yields, in the interaction representation the quantum 
state 



IV'n [t)) ~ exp 



i^i - i^+e-'""* / e'""^C„ (r) z{t)At 



|0), 



(7) 



The approximate expression ([7]) is valid as long as the mean number of photons absorbed per atom remains small. 
With the help of Eqs.®-!!?]) and the functional integral Eq.®, one finds the quantum state of the compound 
atom-cavity system 



I*) 



Vz{t)Vz*{t) 

T 



■ exp 



N 



n=l 



exp 



N 



iJ\z{T)fdT + iJ2'Y^ exp -J J e''^"'^^-^^£{T)£{0)z*{T)z*{0)dTd9 

\o), 



-'^^^n j e'"'^^"^'^Cn{T)z{T)dT exp ia^ J e''^°'-^-*^£{T)z*{T)d 



(8) 



where by \0) we denote the "vacuum state" |0)p JJ^ 10)^^ 

n 

henceforth. Details of the functional integration are pre- 
sented in Appendix A, and the final result Eq. (jA6P is 
given by the following explicit expression 



N 



A{t) 



exp 
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exp 



\o), 



(9) 
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with the coefBcients / given as 
t 

In^ij e'(-o+-")(--*)c„ (r) £(T)dr, (10) 



t T 

In,rn = j j (^) ^(^)^^™ (^) ^(^) 



and A{t) being the time-dependent normahzation. 

Eq. ^ gives the general form of an entangled quan- 
tum state for a quantized cavity field mode coupled to an 
ensemble of N low-excited atoms. The presence of the 
oscillating terms in the integrals of Eqs. pIH) - PT|) guar- 
antee that the parameters /„ and In,m are typically small 
numbers and therefore the probability of excitation per 
atom remains low. However, when uq ~ —w„, in the 
cavity-atom resonant regime, the coefficients /„ can be 
significantly larger as compared to the negligible In^m pa- 
rameters. As we show in Sec lIVl in detail the probability 
for each atom to be in the excited state in this case indeed 
remain small. 

As the next step we establish the connection between 
the parameters in Eqs. pTI)l - (fTT|) . the measured number 
of the cavity photons, and the nilpotcntial's coefficients 
characterizing the entanglement in the atomic ensemble. 

III. METHODS OF ENTANGLEMENT 
CONTROL 

We are now in a position to relate the quantum 
state Eq.® to the nilpotential formulation introduced 
in Ref.Q. Cavity photons can be incorporated into 
the description, with the raising operator af being the 
corresponding 'nilpotent' variable. In particular, for the 
state in Eq.®, the tanglemeter /c, defined as l^*) = 
-^tj^^^ has the form 

JV Af 

/c = a t X! ^nln + ^n^^In,m ■ (12) 

n—1 n,m—l 

According to the entanglement criterion, the first term 
describes the entanglement among the atoms and the 
photons whereby the atoms collectively participate in 
this entangled photons-atoms state via the operator 

N 

X] ^n^n- This term corresponds to the lowest order of 

n=l 

entanglement between the collective state of the atoms 
and the cavity photons. The second term in Eq. p^ is 
concerned exclusively with the entanglement among the 
atoms. In this work we are mostly interested in under- 
standing how the atoms get entangled after the degrees 
of freedom of photons are traced out by a projective mea- 
surement on the "engineered" cavity field state. After the 



measurement of the cavity photon number, the tangleme- 
ter of the multiatomic system undergoes a transformation 
to the generic form 

JV JV 
n,m=l n,m,. ..,; = ! 

(13) 

where higher order terms imply the presence of higher 
order entanglement. For instance, a GHZ state of 
two-levels atoms requires a non-zero term proportional 
to /?^^^ while bipartite entanglement requires only the 
second-order terms to be nonzero. 

In principle, all the (3 coefficients in Ea. (fTn)l are dif- 
ferent and hard to control. We shall therefore consider 
a simpler problem of entanglement between two mul- 
tiatomic ensembles, A and B, each containing N atoms, 
equivalent from the viewpoint of entanglement. 

To describe the entanglement between the two ensem- 
bles, we define for each ensemble a collective nilpotent 

JV JV 

variable: f?^ = ^n^^~^ and = J2 ^r^^^ respec- 

n— 1 n—1 

tively. These variables are not nilpotent in the limit 
N oo since they only vanish in the power N +1. Still, 
the collective operator a\ together with the operators 
(T^ and a% form the su(2) algebra that is a subalgebra 
of the full algebra su(2^) of the A^-atom ensemble. This 
situation belongs to the case of generalized entanglement 
[m , which admits powers higher than two of the creation 
operators. The tanglemeter for two ensembles has then 
the general form 

JV 

fc= Y.f^''A^A)'iKy ■ (14) 

k,i=a 

The presence of higher-order cross terms is associated 
with higher-order entanglement between the two ensem- 
bles, while nonlinear terms that are dependent only on 
one nilpotent variable refer to the entanglement among 
atoms within each ensemble. 

We proceed by listing different methods for entangle- 
ment control, namely, the ways of manipulating the tan- 
glemeter coefficients (3 of Eq.([T3]) for a multiatomic en- 
semble, or that of Eq. (fT4|) for two such ensembles. 

A. Control of the coefficients in Ea. (|12p . 

The P coefficients in the tanglemeter of Eq. (fT3|) di- 
rectly depend on the initial combined atoms-photons 
state in Eq.([9]). Therefore, prior to the field manipula- 
tions and measurement, one can affect the final state by 
controlling the A^(A^-f l)/2 coefficients /„ and In.m in the 
nilpotential Eq.^^. According to Eqs.lUHll-ini), these 
coefficients depend on time-integrals over the prescribed 
time-dependent laser field £{t) and on the coupling pa- 
rameters C'n{t). The latter parameters are determined 
by the cavity geometry, and the individual positions of 
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the atoms inside the cavity. The prospects for reahzing 
this crucial requirement by emerging techniques are 
discussed in Sec. |Vl 

One has complete control over the coefficients In,m, 
when the number of adjustable parameters of the laser 
field and the couphngs exceeds the total number N{N + 
l)/2 of these coefficients. For the simplest control set- 
ting, this implies that a constant field £ is switched on 
and off N{N + l)/2 times, such that each time when the 
field is on, only one of N{N + l)/2 possible pairs of the 
atoms is in the cavity. Given the coefficients /„ and In,m, 
one finds the time intervals when the field is on by solving 
standard linear algebraic equations. This approach holds 
for arbitrary shapes of £{t) and C(t) cast into a super- 
position of N{N + l)/2 linearly-independent functions of 
time, as long as the corresponding linear problem is not 
singular. 

Finally, one can control the relative strength of the 
linear /„ and bilinear In.m coefficients by adjusting the 
relative frequencies, ujq + w„ and loq — LUn, in the os- 
cillating terms of the integrals in Eqs. pO| - (fTT|) . For 
instance, in Sec lIVI control schemes are proposed in the 
cavity-atoms resonant regime, ujq ~ — ct;„, where the bi- 
linear terms In,m in Eg. pip become negligible compared 
to the /„ terms. 



B. Measuring the field in tlie cavity. 

One of the possibilities to control the atomic entan- 
glement is by measuring the number of photons in the 
cavity [l|, [ij, [3 [IBl- The probabilistic outcome after 
having measured d photons yields the non normalized 
state 



1$: 



exp 



■ N 

E 

.n.m— 1 



\0) . (15) 



In the special case of the cavity vacuum, d — 0, the 
nilpotential of the resulting atomic state is already in 

N 



the tanglemeter form /c 



n.m— 1 



The pres- 



ence of only bilinear terms in the tanglemeter indicates 
that the degree of entanglement is not high. Only for 
three d, [l6[ and four two-level atoms 0, [l3| can one 
can prove that this state contain ingenious three-partite 
and four-partite entanglement, respectively. In general, 
for constructing higher-entangled states, one needs either 
to detect d > 1 photons, or to perform certain manip- 
ulations of the field prior to detecting d = photons. 
We dwell on the latter scenario in the next paragraphs 
by considering realistic manipulations on optical cavities. 
The techniques to be presented are hardly practical for 
superconducting microwave cavities where the field of the 
cavity is inaccessible Q. 



C. Displacing the cavity field prior to the 
measurement . 

Let us first apply the field displacement operator D = 
ejl-^'l^e"'*' " e''*"*to the state in Eq.® that experimen- 
tally implies injecting a classical field into the cavity [l^] , 
and then measure the cavity photon number. If d = is 
detected, we keep the projected atomic state, otherwise 
we discard it. This yields the non normalized state 



1$) 



exp 



N 



exp 



■ N 

E 

.n,m— 1 



The local operator exp 



N 



O) . 

(16) 

is nonunitary. In 



combination with the local unitary SU{2) transforma- 
tions it allows one to perform all transformations of the 
group SL{2, C). Therefore, by displacing the cavity field, 
performing local unitary operations and measurements of 
the photon number, we can move the state Eq. (|15p with 
d — 0, along its SL-orbit [l3| thus changing the amount 
of su-entanglement (see Ref.|7'| for more details). The 
operations belonging to the SL{2, C) group can be used 
for entanglement distillation [19|. For three two-level 
atoms (qubits), for example, one can use displacement 
of the field to construct with some probability the GHZ 
state. For more than three atoms the 5i-orbit of the 
state Eq. fTS)) with d = does not contain the GHZ 
state and therefore cannot be obtained by this method - 
the fidelity between the GHZ state and the closest state 
one can construct by this method is decreasing rapidly 
with the number of atoms. 



D. Squeezing the cavity field prior to the 
measurement . 

Squeezing the cavity field state by applying the 
operator exp [(ffS^ — g*&^^ t] prior to the photon num- 
ber measurement offers another possibility to control the 
atomic entanglement. Experimentally, this amounts to 
operating the gas-loaded cavity as a parametric ampli- 
fier. If we start, as earlier, with the state Eq.([S]) 



exp 



afO-l-G 



\0) 



(17) 



with O = E ^n^n and G = ^ a^a^^In,m and de- 

n—l n,m — 1 

tect zero photons in the cavity after the squeezing, the 
reduced state of the multiatom ensemble reads 



*"'')*e^to+G|o)p|0) 



*"'')'e^tOe<5|0)p|O) 



^5(3+6)^t~g*3+^t+G 



|0),|O) 



(18) 



For a real squeezing parameter g, the non normalized 
state Eq.lIH]) can be rewritten as |<I>) = exp(CO^-f G) 1 0) , 
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(see Appendix IB] for details). The presence of the square 
of the operator O in the atoms' nilpotential impHes that 
the atoms can be entangled with each other even if zero 
photons are detected and the cavity-atoms frequencies 
are tuned to the resonant regime where G is negligible. 
In particular, this can be done for two ensembles as shown 
in Sectdl 



E. Creation of highly entangled atomic states by 
field nonlinearity 



Apart from displacement and squeezing there exists 
another tool for the field manipulation: a Kerr-nonlinear 
gas (with appreciable nonlinearity at the cavity-mode fre- 
quency Wo) can be introduced in the cavity after the 
atoms have passed through it. This type of coupling 
results in a nonlinear dependence of the cavity energy 
on the number of the cavity photons. The presence of 
the medium can also couple an external laser field at a 
frequency lol with the cavity field, thus inducing a mul- 
tiphoton cavity excitation. 

To be more specific, let us consider a symmetric Kerr 
medium with the nonlinear polarization P = ^E'^, where 
the electric field E consists of the classical laser field 
of large amplitude £ and the quantum cavity field ^ 
(a f +a) . In the rotating frame defined by the unitary 
transformation Uq — exp [— io^^a f at] , the interaction 
energy PE integrated over the cavity volume V yields 
the Hamiltonian 



{LUo-ujL)a'fa + K (afa) +(at+a)f^ , (19) 



with K ~ ^V. By adjusting the laser frequency lol and 
amplitude £, one finds a multiphoton resonance with the 
cavity photons and creates superpositions of the cavity 
photon states that differ by a large fixed number of pho- 
tons. If one detects the cavity in the vacuum state af- 
ter the initial state Eq.® has been subject to such a 
transformation, the atomic ensemble turns out to be in 
a highly entangled state. This can become a GHZ-state 
for certain values of the parameters, as it will be shown 
in the next section (Sec. IIVB|) . 



IV. APPLICATIONS 



Having presented the main ideas of the control schemes 
we illustrate the suggested techniques for specific states, 
such as the A^-atom GHZ and the A^-atom Dicke state 
with M excitations. The W state of A^-atom is a Dicke 



state with M = 1. The state vectors of these states read 
I00...0) + 



\GHZ) 



V2 



\W) = \1,N)^ = 



A^ 



M 



(20a) 



respectively, where {Pk} denotes the set of all distinct 
permutations of atoms. The corresponding tanglemeters 
are of the form 

N 

fGHZ=^cj+o+ ...a+- fw^ E op}- (20b) 

M JV M N 

j = l i=M+l k,m=l i.j=M+l 



M N 
1=1 j,k,...l=M+l^ 



M 



These states are often discussed in the context of quan- 
tum information processing. Moreover, a prominent ex- 
perimental realization of quantum continuous variables is 
obtained, when the quantum states of an asymptotically 
large ensemble are treated collectively, admitting opera- 
tors with asymptotically continuous spectra (20| . There- 
fore one of our examples is concerned with the engineer- 
ing of entanglement between two atomic ensembles. 

We note that all the constructions to be presented are 
probabilistic in nature and the desired state is attained 
only when the cavity field is detected in a specific Fock 
state. Yet, the success rate can be appreciable. Even 
though the methods presented so far concern the general 
case where all the / coefficients Eas. p^ - ([TT|) are nonzero, 
we shall next consider applications in the resonant cavity- 
atoms regime where the In,m are negligible, while the /„ 
are not. 



A. Constructing a Dicke state with M excitations 

Consider A^ identical atoms that are sent through the 
cavity of resonant frequency wq ~ — w„. If all the atoms 
are manipulated equivalently, the combined wave func- 
tion Eq.® yields the state 



1 

1 



exp 



N 



n=l 



(21) 



with identical field-atom coupling coefficients It = c, and 
vanishing pairwise coefficients = 0. The normaliza- 
tion factor A is 



N 



iim 



(22) 
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while the excitation probability for each atom is 1/A''. 
Therefore, for large N, Eg. ipT]) is consistent with our 
initial assumption of low excitation per atom. 

When a measurement of the cavity photon number is 
performed, the result d — M is obtained with probability 



P(M, c) 



M\N\ |„|2M 
{N-Aiy. 1^1 

(JV-i)! |C| 

i=l 



(23) 



leaving the atomic ensemble in the \M,N)jy state. Ex- 
perimentally one can adjust the coupling coefRcient li — 
c such as the probability of attaining the desired Dickc 
state becomes maximum. In Fig. 1 we show for the 
ensembles of = 10 and = 19 atoms, the proba- 
bihty P{M,c) Eq.dlJ) for different values of the cou- 
pling coefRcient |c| and for different excitations M. We 
note that the Dicke states \M, N) jy and |iV - M, 7V)^ al- 
though they are equivalent up to local operation on each 
atom, the maximum probability for the creation of each 
of them might be considerably different. Therefore, one 
should choose to construct the one that acquires highest 
probability P. 





FIG. 1: (Color online) The probability P{M,c) for obtaining 
the Dicke state \M, N)^^ as a function of the amplitude of the 
coupling coefficient |c|. We consider two ensembles consisting 
of A'^ = 10 and A'^ = 19 atoms, and a different number of 
excitations M. 



Construction of the GHZ state of an atomic 
ensemble 



We begin with the observation that the state Eq. ([2T|) . 
corresponding to the lowest possible order of entangle- 
ment between the individual atoms and the cavity field, 
may still yield a highly entangled state of the atomic en- 
semble after being projected onto a linear combination of 
the field vacuum and a highly excited Fock state. Let us 
take the linear-combination field state 



\F)^B\0)+C\N) 



p ' 



(24) 



where A" is the total number of atoms. Then this projec- 
tion can be shown to result in the GHZ state. In fact, by 
casting Eg. pT]) in the Taylor series 



1 V^ (Qt) sr-^+T 



\o) 



(25) 



= 1 \n=l 

one immediately finds the projection 

N 



{F\^) 




\0) 



For 



N 



B* 



(26a) 



(26b) 



this indeed yields the GHZ state with the tanglemeter 
fcHZ of Eq.(l20b|. 

Detection of the field state cannot however be 
directly performed as a probabilistic projection on the 
photon number basis and requires a manipulation of the 
field prior to such a measurement. One way to perform 
such a manipulation is to expose the system to the Kerr- 
nonlinear interaction Eq. l|19p with the parameters chosen 
such that only the corresponding states with and A^ 
cavity photons are resonant and coupled. This condition 
reads 



0, 



(27) 



while the corresponding matrix element of the multi- 
photon transition between the photon vacuum and A^-th 
Fock state has the form 



Vr 



ON 



73N 



N 

n (nn^ 

n=l 



ujLn) 



N 

/Mn {n-u;L/K) 

n=l 



(28) 



The interaction Eq. (|19p must be on during a time 
given by the condition 



tan {tKerrVoN) 



C* 

'b* 



N 



n In 



(29) 



7 



Then, if the cavity field is detected in the vacuum state, 
the projection onto the state \F) is performed. This re- 
sults in the GHZ multiatomic state. 



C. Entangling two ensembles 

We now consider entanglement between two ensem- 
bles of atoms, A and B, each consisting of N identical 
atoms. Each ensemble is treated as a single element and 
we are exclusively concerned with the collective entan- 
glement between these multiatomic elements expressed 

N , 

in terms of the collective operators — and 

n=l 

N 

~ ^ We guide the first ensemble of atoms 

through the cavity and then repeat the procedure for 
the second ensemble. If the excitation per atom remains 
small and the cavity is tuned on resonance with the atoms 
frequency, then according to Eq.® the combined state of 
the two multiatomic ensembles and the cavity field reads 

1^-) = le'^st(s++s+) ^3Q) 

Entanglement between the two ensembles is implied by 
the presence of cross terms ((?^) (o^b) nilpoten- 
tial, can appear as a result of the field manipulation prior 
to the measurement of the cavity photon number. In 
accordance with the results of Sect lllll the lowest-order 
cross term with k = I = 1 can be generated from the state 
Eq. ([5(I)l . by squeezing the cavity field followed by the de- 
tection of the cavity vacuum. More specifically, after the 
squeezing operator S — exp (gS^ — 5Q?^^) has been ap- 
plied to the state Eg. ([50]) for time t, and zero photons 
are detected, the ensemble state adopts the form 

|$) = -^exp(C02+ 77)10) , (31a) 

where 

r = —===, ri = gt, C = 2t(inh[gt]- gt , (31b) 
Vl + e^^ 

and the value of the parameter fj,, given by Eq. (fTO|) . is 
identical for all the atoms. 



V. DISCUSSION 

We have presented a powerful formalism for the ana- 
lytical description of multiatomic ensembles interacting 
with quantized cavity fields for the case where the atoms 
have low excitation probability and the effect of decoher- 
ence is ignored. By combining the analytical functional 



integration with the multipartite entanglement descrip- 
tion via nilpotent polynomials we have been able to pro- 
pose general schemes for the probabilistic control of en- 
tanglement among multiple neutral atoms. These tech- 
niques consist of manipulations on the cavity field fol- 
lowed by a projective measurement on the photon num- 
ber state and therefore are experimentally feasible only 
in optical cavities. We presented several applications in 
the regime where the frequency of the cavity mode is 
resonant with the atomic transition frequency. 

The proposed control scheme presumes the resolution 
of various experimental problems posed by state prepa- 
ration and detection, most importantly the determina- 
tion of atom number N in the cavity. There is continu- 
ing progress towards the measurement of trapped atom 
numbers, suggesting that this goal may be achieved be- 
fore long. The following procedure may be conceived of: 
(?) slowly trapping ground-state atoms in the cavity (for 
which several options exist, including the adiabatic con- 
veyor - belt technique (**) counting the trapped 
atoms by their resonance fluorescence or other optical 
techniques; (Hi) impressing Stark or Zeeman shifts to 
make atoms at different positions spectrally distinguish- 
able and thus addressable by the control Raman cou- 
pling; (iv) switching (on and off) the Raman control field 
in Hamiltonian ^ for selected (spectrally addressable ) 
atoms for the required interaction time. This procedure, 
although challenging, may in the long run allow field - 
atom entangling manipulations without actually sending 
atoms in and out of the cavity, but rather making them 
disappear and reappear by purely optical manipulations. 

Beyond the ability to realize a broader variety of en- 
tangled multiatom states than existing states, the main 
merit of the present analysis is that it has been exhib- 
ited within the unified framework of the general nilpotent 
formalism. 
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APPENDIX A: FUNCTIONAL INTEGRATION 

The straightforward way to calculate the functional in- 
tegral of Eq.® is by transforming it into a Gaussian 
functional integral. This requires first to find the sta- 
tionary solutions z*{t), Zs{t) that correspond to the ex- 
tremum of the functional exponent, i.e., the action 5, 
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t t 



and then to evaluate the integral for new functional vari- 
ables {z* (t) + z*{t), z (r) + z,{t)}. 



S = -i f\zir)\'dr + iY^'^t 

"=i 

N * 

f e'-"(--*)C„(r)z(T)dr 

n=l { 
t 

+ iSt / e'"«("~*)f(T)z*(r)dT 

e''^«(^-^)£:(r)f(6l)z*(r)z*(6»)drd6', (Al) 



T 



By variation of the action S 



O^SS= / dr 



N 



-iz*(r)~i^a+C„(r)e-'-(*-^) 



iz(T) -iate'"°'^-*)f(T) +£:(t) y e'"°(^-^)f(6')z*(6i)d6' + f(T) j £(6i)z*(6i)e''^°(^-^M6' 

T 

one obtains the stationary solutions 

N 

<W = -E^nC'„(r)e-'-"(*--' 



(5z*(t) (A2) 



n=l 



N I 



e)+ic.o(9-r)^^ . (A3) 



r 



The substitution of the new variables into Eq. ljAip . 
separates the action 5* into two parts 



S = S{z (r) , z* (r)) + 5(z, (r) , zl (r)), (A4) 
a "quantum" part and a "classical" one. The quantum 



contribution of the action does not contain quantum op- 
erators in our case and the functional integration over 
this gives just a phase that can be ignored. What is left 
then is to evaluate the classical contribution to the func- 
tional integral by substituting Sd — S{zg (t) , z* (r)) 



N N 'j, 

= i E - ^ ^ y e'(-+-)(-*)C„ (r) £(r)dT 



n=l 
N 



■ E ^n^m / / f(r)£(0)C„(0)C„, (r)e-'("("«-""'+*(""+"'")-^(""+"»»dMT. (A5) 

n,m— 1 
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into Eq.®. The final expression for the wave function describing the time-evolution of the combined cavity-atoms 
state is 



Ait) 



exp 



exp 



N 



N N 'j, 

i^^t + iSt^a+ / e'(-+-)(-*)C„(r)f(r)d7 

„ — 1 ^ — 1 



n—1 n—1 

t T 



n,m— 1 







\0) . (A6) 



where A{t) is the normalization factor whose derivation 
is presented in the following. 



1. Normalization factor 



By definition the square of the normalization factor, using the notation of Eqs.((9|)- (fTT|) . is 



|^(t)r = (0|exp 



JV 



iV 



exp 



N 



N 



\0) , (A7) 



or, equivalently, after expanding the terms involving the operators a and af , 



|A(t)|2 = (0|exp 



N 



T* 



exp 



N 



N 



n—1 m— 1 



exp 



AT 



_n,m — 1 



|0) 



(A8) 



Under the assumption of low excitation probability per atom, the collective operators approximately commute and 
the Eq. ljASp may be rewritten in the form 



\A{t)\' ^ {0\e^V 



N 



N 



N 



N 



n—1 m— 1 



|0) 



The expression (jA9p can be written in a simpler way as 

\A{t)\' = {0\exp[yMy]\0) , 
by introducing a Hermitian matrix M and a complex vector y 



M 



J* 



N 



\InI*n I 



'Af,l 2 



N,l 



h. 



N 



\ 



Il,N 

\InI*n I 



•n 



(A9) 



(AlO) 



(All) 



At this point, we can employ the identity 



dx dx exp [— {x + w) A{x + w)] — 



(2^) 
det^ 



(A12) 



that holds for complex vectors x and w, to prove that 



det A 



(2^) 



N 



dx dx exp [— — zx — xz] 



(A13) 
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for A, a Hermitian matrix and z, a vector defined as z = Aw. 

I 

With the help of the identity (|A13|) we can rewrite Eg. ljAlOp as 

|^(i)|' = (0|exp[fMy]|0) 
det Af-i 



2^ . dx^^dr^^ {0\ exp [-xM~^x -yx- xy] \ O) 
dx^^dr^^ exp [-xM~'^x + [x + x)V{x + x) 



dx'^^dx^^exp 



( x x) B 



(2^) 
det 

_ det M^i 

_ det M-i 
~ detB ' 

where we have assumed that Xi small, an assumption valid in the chosen regime, and 

r 



V V 

V V-M-^ 



(A14) 



(A15) 



APPENDIX B: SQUEEZING OPERATORS 

Having projected the squeezed state onto the vacuum 

1$) = (OpI c(^("+°)'-^*"'')* |0p) e^ I O) (Bl) 

we would like to explicitly calculate the state of the atoms 
in the special case where g = g* ■ If the position- 
momentum representation of the ladder operations is em- 
ployed, 

1 5 , 



V2 dx' 



(B2) 



the problem is reduced in solving the Shrodinger equation 
for the Hamiltonian 



H ^i\g[a + 



ga 



t2 



{y2g^ + gx^^+^(^/2g^x + g + g^^'') (B3) 



and with initial condition 



(x|Op) 



(B4) 



If g is a real number, the differential equation to be 
solved is of first order both in position and time and thus 
can be easily handled. The final result is 



I $) = r exp 



exp 



G 



\0) 



(B5) 



CO' +?7 

where t is time period that the squeezing operator is ap- 
plied and the parameters are 



r = V27r/Vl + e25*, rj = gt, ( = 2 tanh[gi] - yt. (B6) 
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